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A theory of the proximity effect in normal metal - multiband superconductor hybrid structures is formulated 
within the quasiclassical Green's function formalism. It is shown that the existence of multiple superconducting 
bands manifests itself as the occurence of additional peaks in the density of states in the structure. The interplay 
between the proximity effect and the interband coupling influences the magnitudes of the gaps in a supercon- 
ductor in a non-trivial way. The developed theory is applied to the calculation of supercurrent in multiband 
Superconductor - Normal metal - Superconductor Josephson junctions with low-transparent interfaces, and the 
results are compared with the predictions for multiband tunnel junctions. 
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The proximity effect is the phenomenon that a supercon- 
ducting order parameter can penetrate from a superconductor 
(S) into a normal metal (N), or another superconductor (S') 
with a critical temperature T c s< < T c s, over a distance of the 
order of the coherence length, inducing a minigap in N or S. 
This phenomenon is well understood, both in terms of An- 
dreev reflections as well as in terms of microscopic Green's 
functions, 1 i 2 , 3 ,ViV, 8 . 9 

It is not known, however, how the proximity effect will 
manifest itself when multiple pairing potentials are present 
in the superconductor. This question has become relevant 
now that multiband superconductors are coming into practical 
use. The most clear example of a multiband superconductor 
is MgB2, for which the experimental and theoretical evidence 
for the coexistence of two gaps is overwhelming. 10 The multi- 
band nature of the superconductivity in MgB2 is theoretically 
well explained 11 by the qualitative difference between differ- 
ent sheets of the Fermi surface, together with the large dis- 
parity of the electron-phonon interaction. Therefore, in this 
paper, the question is addressed how the multiband nature in- 
fluences the proximity effect. E.g. what will the density of 
states look like in a SN bilayer, where S is a two-band super- 
conductor? 

Josephson and quasiparticle tunneling in hybrid structures 
containing multiband superconductors have been investigated 
theoretically in Ref. [l2j and applied to the calculation of the 
total Josephson current in a SIS two-band Josephson tunnel 
junction. For all-MgB2 devices, high-quality tunnel barriers 
are not available yet, and realizing SNS structures is an at- 
tractive alternative, of which first systems have been realized 
already. 13 In this paper, the theory of the multiband proximity 
effect is applied to the calculation of Josephson current in SNS 
structures having two-band S electrodes. The practically inter- 
esting SINIS case is considered, where a non-ideal interface 
transparency is taken into account. Predictions are made for 
Josephson devices based on MgB2 and compared with those 
for MgB2-based tunnel junctions. 

In this paper, we will use the quasiclassical Green's func- 
tion formalism in order to describe electrical transport in SS' 



hybrid structures, where S' is a single-band superconductor 
while S is a multiband superconductor. We will restrict our- 
selves to the limit of diffusive transport, which is justified if 
Is.S' <^ £s,S'> where l$,s' an d £,s,S' are the electric mean free 
path and coherence length of the S and S' materials respec- 
tively. In the dirty limit, the Green's functions in the S' metal 
are given by the standard Usadel equations. 14 In the S metal 
in the regime of vanishing interband scattering, as is the case 
for MgB2, 15 the Usadel equations take the following form 16 

^[(G§) 2 $^]'-*§ = -A ai (1) 

A a =2nT V A a/3 ^^. (2) 

* — ' U! 

P,LJ>0 

Here, a and j3 are the band indices, e.g. a, f3 = 1, 2 in the 
two-band case (later we will use the band indices a and 7r for 
MgB2 specifically), A a is the pair potential, Gg and <f>g are 
Green's functions, 5 u> — irT(2n + 1) are Matsubara frequen- 
cies, Dg is the diffusion coefficient, and K a p is the matrix of 
effective coupling constants. The prime denotes a derivative 
with respect to the coordinate x in the direction perpendicular 
to the S-S' interface. 

Equations Q and (|2ji must in general be supplemented by 
boundary conditions. Zaitsev 17 derived boundary conditions 
to the quasiclassical Eilenberger equations at the S-S' bound- 
aries in the clean limit, which were further simplified in Ref. 5 
in the dirty limit. These boundary conditions have to be mod- 
ified when S is a multiband superconductor. 

In the limit of small interband scattering a multiband super- 
conductor may be represented by separate groups of super- 
conducting electrons which interact with each other only indi- 
rectly, via selfconsistent pair potentials in the bulk. Therefore, 
for the derivation of the boundary conditions for the Usadel 
equations, one can apply a similar procedure to that used in 
Ref.|5|in the single-band case. In the multiband case, the set 
of interface parameters, 7° and 7^, describing the proximity 
effect, should be introduced for each of the bands. 



2 



The first boundary condition relates the current from the S' 
metal side at the S-S' interface, ctG 2 $', to that from the S 
side, J2 a o-s(Gg) 2 ($g)'. Therefore we have 
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Here, a = 1/p and erg = l/pg are the conductivities of the 
S' layer and the respective bands of the S metal, D is the dif- 
fusion constant in S' and T c $ is the critical temperature of S. 
The ratio between the parameters 7" for the different bands is 
mainly determined by the relation between the diffusion con- 
stants, Dg. In the case of MgB2, the 7r-band is generally con- 
sidered to be more dirty than the a-band, 15 i.e. Dg <§C Dg. 

The second boundary condition relates the gradient of the 
Green's function $ near the S-S' interface to its jump at the 
interface due to the finite interface resistance, 



= E § (*> 

IB 



(5) 



where 7^ = R B /pC R% are the components of the specific 
interface resistance, describing the tunneling of an electron 
across the interface into the corresponding conduction band. 

In order to obtain the resistances R B , we have to evalu- 
ate the effective junction transparency components. It was 
first pointed out by Mazin, 18 that the normal state conductance 
i?" 1 , in the limit of a specular barrier with small transparency, 
is proportional to the Fermi-surface average ( x Nv 2 ') a , where 
N is the density of states and v the Fermi velocity. In Ref.[l2l 
it was further shown that the normal state resistance compo- 
nent of tunneling into band a of S is given by the contribution 
of the electrons in band a to the squared plasma frequency 

(ojp) , which can be obtained from first principle calcula- 
tions. For MgB 2 , the ratio R B /R B = i^p/^pf is 2 and 
100 for tunneling in the direction of the a — b plane and c-axis 
respectively 

In the case of a SS' bilayer, the Usadel equation ([0 needs 
to be solved in the S as well as in the S' layer, together 
with the self-consistent determination of the pair potentials 
in S and S', Eq. A general numerical method, using 0- 
parametrization, $ = wtan0 and G = cosG, is described for 
the single-band case in Ref.jg. Here, we extended this method 
by applying the new boundary conditions, Eqs. (0 and (|5j- 
The density of states at energy E can be obtained by applying 
an analytical continuation, u = —iE, to the Usadel equa- 
tions and the boundary conditions and solving the numerical 
scheme in the complex energy plane. 

The numerically obtained dependence of the pair-potential 
on position is presented in Fig. [2 for the example in which 
the coupling constants are taken as calculated for MgB 2 in 
Ref. llSJ The parameter values are indicated in the caption. 
For temperatures above T c g> (solid lines in Fig. [0, it can 
be seen that the pair-potential in S' increases towards the in- 
terface, while A CT decreases, as expected in analogy with the 
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FIG. 1: Pair-potential as function of position for a SS' bilayer at 
T = 0.5T c s. The parameters of the bilayer are r y"' n = 1, jg = 2, 
7S = 1, ds/^s = ds'/fs' — 10, and the coupling constants in 
the S-layer are chosen as expected 19 for MgE^: An = 0.81, A22 = 
0.278, A12 = 0.115, and A 2 i = 0.091. 



proximity effect in the single-band case. The decrease in A CT 
towards the interface can be explained by the relatively strong 
coupling between the a and it bands. By decreasing the inter- 
band coupling constants and by decreasing the interface sup- 
pression parameters, one can obtain the opposite regime, in 
which increases towards the interface. For relatively large 
values of T c s>, and for temperatures below T c s', we even pre- 
dict an increase in A^ towards the interface, as illustrated by 
the dashed line in Fig. ^ This effect is further enhanced when 
the cr-band is decoupled from the S' -layer, which is the case 
for example when the interface normal is parallel to the crys- 
tallographic c-axis of MgB2, due to the vanishingly small ratio 
R B /R B in that case. 

As an example, in Fig. |2]the results of a calculation of the 
density of states in a SN-bilayer are presented. In the con- 
sidered case, the bulk energy gaps in a two-band supercon- 
ductor are not too different. As is seen from the figure, the 
density of states in the N-layer has three peaks: the lowest en- 
ergy peak corresponds to the proximity induced minigap and 
the two other peaks correspond to the bulk energy gaps in the 
two-band superconductor S. The existence of a minigap is a 
characteristic feature of the proximity effect in a SN-bilayer 
in the dirty limit, as was studied in detail in the single-band 
case in Ref. 6. As we can see, the minigap persists in the 
two-band case as well and its magnitude depends on the pa- 
rameters of the interface, thicknesses of the N and S and the 
values of the bulk gaps in the superconductor. 

The next step in investigating the influence of multiband 
superconductivity on the proximity effect, is to study super- 
currents in multiband proximized structures. We will consider 
double-barrier structures consisting of two S electrodes cou- 
pled by a normal metal N. As a model system we use a SI- 
MS double-barrier hybrid structure, since in practical devices 
interface potential barriers are always present at the S-N in- 
terfaces, either originating from a Fermi-velocity mismatch, 
degradation of surface layers, or artificially deposited oxide 
barriers. 
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FIG. 2: Normalized density of states in a proximized SN-bilayer at 
several positions in the bilayer ( 1 — 4, as indicated in the inset), where 
S is a two-band superconductor. The parameters of the bilayer are 
7 1 ' 2 = 0.1, 7jj' 2 = 5, ds/^s = 10, d/£ = 1 and the coupling 



constants in the S-layer are chosen as An 
A21 = 0.1. 



0.5, A 22 = 0.4, A12 



If the conditions of the dirty limit (electron mean free path 
^ < are fulfilled in the N interlayer, than the station- 
ary Josephson effect in the structure can be analyzed in the 
framework of the Usadel equations by the method developed 
in Refs. I5landfl7lfor the single-band case. We assume that the 
interface transparencies are small enough such that the condi- 
tion 1 + 7§j 2 3> 7" 2 holds at both NS interfaces (here and 
below we drop the subscript S). In this case, the suppression 
of superconductivity in the S layers is weak and the Green's 
functions in the electrodes near the interfaces, G" 2 an d $121 
are equal to their bulk values. To calculate the supercurrent, it 
is sufficient to consider Eq. at the two interfaces, giving 
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For simplicity, we will consider symmetric junctions where 
G" 2 = Gg and Jbi 2 = 7b> an< ^ where the functions $" 2 
are related to the phase shift ip across the junction by 2 = 
A a exp (±iip/2). Further, we consider purely normal N-layer 
with vanishing pair potential A = and restrict ourselves to 
considering the limit of a small interlayer thickness, d -C £. 

In the limit d <C £, there are two characteristic frequencies 
ill. 2 in the Usadel equations dll2> . At u> < fii = t:T c s^ ^ 
il 2 = irT c s we can neglect all nongradient terms in the Usadel 
equation. Hence, [G 2 ^'] = 0, and in the zero approximation 
on d/£ one obtaines that all $ functions are spatially indepen- 
dent constants, $ = A. In the next approximation we have 



* = A + B 1 + A-£- t 



P 2 = 



nT cS G 
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From the boundary conditions and by taking into account that 



FIG. 3: Normalized density of states in the interlayer of a SINIS 
double-barrier structure, where S is the two-band superconductor 
MgB2, and the phase difference over the junction is n/2. The den- 
sity of states is shown for small 7bm (dashed line: 73 M = 0.2, 
Ibm = 0.1) and large 7s m (solid line: ^ BM = 20, y BM = 10). 
For comparison, the inset shows the density of states in the inter- 
layer of a SINIS junction with single-band superconductors (solid 
line: 7bm = 2, dashed line jbm = 2.10 -3 ). 



in our model 2 = A a exp (±i</j/2), we finally will have 
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where Jbm = Ib^/C V 2 — cos 2 ((p/2), and 
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The density of states N(E) — Re(G) in the interlayer of 
the double-barrier junction can now be found from an analyt- 
ical continuation of Eq. (|9jl to real energies, u> = —iE. The 
results for the two-band case are plotted in Fig. [5] The known 
density of states for a single -band SINIS junction 21 is shown 
in the inset. For •jbm <C 1, the single-band results show a 
peak in the density of states at Acos(</>/2), while the density 
of states in the two-band junction in this regime is predicted 
to have a peak at a value that is even lower than A,rCOs(0/2). 
For larger values of 'Jbm, the density of states shows three 
peaks: at the minigap and at A T and A CT , in analogy with 
the two peaks in the density of states of a single-band SINIS 
junction. 

Substituting Eq. into the supercurrent expression, 



I = cr27rTIni V -^G 2 <T$', 



(ID 
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and taking $*$' in lowest order equal to A*B, we obtain 
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dA 2 sin(vj) 
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FIG. 4: I c R h 



regime of 7^ 



for double-barrier MgB 2 SINIS junctions in the 
■C 1 (solid lines), compared to I c Rn for MgB 2 
SIS tunnel junctions' 2 (dashed lines). The total I c Rn ofa — b plane 
MgB2 junctions is an average over all bands, while c-axis junctions 
only contain a 7r-band contribution. 



A generalization to take boundary asymmetry and a finite A 
in the interlayer into account can be made straightforwardly. 

In the two-band case in the limit 7^ M — > 00, which is for 
example the case for tunneling in the MgB2 c-axis direction, 
the normal metal is only proximized by the 7r-gap of electrode 
S and Eq. (I12t gives 



I = 



ttT 
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which has been previously obtained 5 ^ for SINIS junctions 
with single -band superconductivity in S. It also follows from 
Eqs. H2\ and (II 31 that the critical current of a single-band 
SINIS junction always exceeds the critical current of a two- 
band SINIS junction with a vanishing second gap, = 0. 



The temperature dependence of the critical current can now 
be calculated for SINIS Josephson structures for different ori- 
entations of the crystallographical axis with respect to the in- 
terface normal. The gap functions, A 7r . <T (T) and the ratio 
1b/1b' follow from band structure calculations. 12 The re- 
sults are shown in Fig. |4]for vanishingly small 73 m, and 
compared to the calculation results for SIS junctions.— The 
full specific interface resistance of a SINIS junction, Rjy = 
RbRb/[Rb+Rb]- 14 is clearly seen that the critical current 
of SINIS junctions is larger than in SIS structures, practically 
in the whole temperature region, as is the case for single band 
superconductors. 5 ' 20 At low temperatures the I c Rn product 
can be as large as 5.2 mV when only the 7r-band contributes 
to the current and close to 7.3 mV when the sum over dif- 
ferent band contributions can be taken into account, as is the 
case for tunneling in the direction of the a — b plane. The nega- 
tive curvature of I c Rn(T) is a direct consequence of the two- 
band nature of superconductivity and is absent in I c Rn(T) of 
single-band SINIS junctions in the regime of small jbm- 20 

In summary, we have formulated a microscopic theory of 
the proximity effect in hybrid structures based on multiband 
superconductors in the diffusive limit. We have shown that the 
existence of multiple superconducting bands manifests itself 
in the proximity effect between a normal metal and a super- 
conductor as the occurence of additional peaks in the density 
of states at the normal metal side. The interplay between the 
proximity effect and interband coupling determines the gap 
magnitudes at the interfaces. The supercurrent in multiband 
proximized Josephson junctions was calculated and compared 
to known single-band results and predictions for multiband 
tunnel junctions. 
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